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Motivation 


• Shock-Wave / Turbulent Boundary Layer Interactions (SWTBLIs) are 
pervasive to all supersonic vehicles - external aerodynamics, inlets, 
isolators, etc. 

• Accurate prediction of SWTBLIs with CFD remains elusive. 

• Results from the 2010 AIAA SWTBLI workshop indicated deficiencies in 
RANS turbulence models for shock-separated flows. 

e Some promise using LES/DNS, but these methods are not yet ready for 
engineering applications. 
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UFAST SWBLI Test Case 

2010 A! AA SWTBLI Workshop 

• Mach 2.25 flow approaching SWTBLI region. 

• Several RANS and LES (including hybrid RANS-LES) solutions submitted; 
most widely used RANS turbulence models were SST and SA. 
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Motivation 


• 2010 AIAA SWTBLI workshop results: Most popular turbulence models 
utilized were MenterSST and Spalart-Allmaras. 

• SST and BSL produce results on either side of experimental data for this 
case (UFAST) and most other SWTBLI cases we have examined. 

U FAST Mach 2.25 SWTBLI Test Case: 


Experiment 
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Spalart-Allmaras 
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Presentation Outline 


Motivation 

Comparison of BSL and SST Models 
Investigation of Shear Stress Limiter 
Incompressible / Low-Speed Test Cases*: 

1. Flat Plate - zero pressure gradient boundary layer. 

2. Driver Axisymmetric Diffuser - adverse pressure gradient with separation. 

3. Backward Facing Step 

• SWTBLI Cases: 

1. UFAST- Mach 2.25 

2. Schulein - Mach 5 

3. HIFiRE Scramjet Flow Paths - Mach 5.8 and 8.0 flight conditions 

8 Conclusions 


* Grids and boundary conditions for incompressible test cases are taken from the Turbulence 
Model Benchmarking Working Group (TMBWG) website: turbmodels.larc.nasa.gov. 


Menter SST and BSL Models 

• Baseline (BSL) model combined: 

1. Wilcox 1988 k-co model - inner model. 

2. Jones Launder k-s model - transformed to k-co equation form - outer model. 

- F x function transitions from inner model to outer model at approximately 2/3 

boundary layer thickness. 

• Shear Stress Transport (SST) model is an extension of BSL model: 

1. Diffusion coefficient, a k , changed from 0.5 to 0.85 for inner model. 

2. Limiter placed on turbulent shear stress to not exceed 0.31 x turbulent kinetic 
energy (TKE). The motivation is to account for "transport of shear stress." 

• The limiter originates from the observation of Bradshaw and others 
that -u V does not exceed 30% x TKE 

• "Townsend structure parameter" = a 1 = -u V / k 

• SST model setting a x = 0.31 works well for low speed cases including 
mild adverse pressure gradient flows. ...optimal for higher speed? 
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Menter Two-Equation k-co Baseline (BSL) Model 
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Menter Two-Equation k-co “Shear Stress Transport” 

(SST) Model 
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Investigation of Shear Stress Limiter 

• The original SST model has become one of the most widely RANS used 
turbulence models. 

• Shear stress limiter sets a 1 = 0.31 = -uV/k, using observations of 
Bradshaw, Townsend, and others -- for zero pressure gradient and mild 
adverse gradient flows. 

® This limiter is only active in inner % of boundary layer - via F 2 function. 

• Experimental data (UFAST, Smits et al, others) shows -uV/k exceeds 0.31 
in SWTBLI flows. 

• Others (Wilcox, Tan and Jin, Edwards) have investigated values for a 1 
greater than 0.31. 

• This work: 

1. Investigate a range of a 1 from 0.31 to 0.40. (larger values for a x are very 
similar to BSL) 

2. Examine details of experimental turbulent measurements alongside 
computations to determine appropriate value(s) for a 1 in SWTBLI flows. 
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Mach 0.2 Boundary Layer 

Grid, BC's from TMBWG website. 


Wall Skin Friction: 
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Zero Pressure Gradient 


Mean velocity (Re x = 4.3e6): 
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Mach 0.2 Boundary Layer 


Eddy Viscosity (Re x = 4.2e6): 



Zero Pressure Gradient 


Shear Stress (Re x = 4.2e6): 
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Driver Axisymmetric Diffuser 

• Grid, BC's from TMBWG website. 

• Original flow case was axisymmetric diffusing geometry within rectangular 
wind tunnel. Menter (AIAAJ. 1994) defined an axisymmetric streamline to 
simplify calculations. 

Wall Skin Friction: Wall Static Pressure: 
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Driver Diffuser - Turbulent Shear Stress 



Experiment: 


SST: 
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Driver Diffuser - Structure Parameter = -u'v'/k 


Experiment: 


SST: 


SST, a ± = 0.34 : 


SST, a 1 = 0.355: 


SST, a x = 0.37: 


BSL: 
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Driver Diffuser - Turbulent Kinetic Energy 



Experiment: 


SST: 


SST, a 2 = 0.34: 


SST, a 1 = 0.355: 


SST, a x = 0.37: 


BSL: 
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UFAST SWBLI Test Case 

2010 AIAA SWTBLI Workshop 

• Mach 2.25 flow approaching SWTBLI region. 

• Several RANS and LES (including hybrid RANS-LES) solutions submitted; 
most widely used RANS turbulence models were SST and SA. 
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Structure Parameter - UFAST Experimental Data 

• Turbulent shear stress exceeds 0.35 x TKE at beginning of interaction 
region and in region where boundary layer reattaches. 


Mean velocity: 



Structure parameter: 


20 h 



-u'v'/k: 0 0.1 0.2 0.3 0.4 
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Structure Parameter - UFAST Experimental Data 
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UFAST Velocity Contours 


Experiment 



MenterSST k-co 
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UFAST Velocity Contours 



Experiment 


MenterSST k-co, a x = 0.34 


20 


■ ■ WM 

20 



U (m/s): -50 50 150 250 350 450 550 



U (m/s): -50 50 150 250 350 450 550 



'280 


300 


320 x (mm) 340 

Menter SST k-co 


360 



'280 


300 


320 x (mm) 340 

Menter SST k-co, a 1 = 0.355 


360 



Menter BSL k-co 


Menter SST k-co, a 1 = 0.37 



UFAST Velocity Profiles 
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UFAST Velocity Profiles 



X = 320 mm 


X = 340 mm 
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Mach 5 SWTBLI - Schulein 


Only wall shear stress data available for this case. 




Hypersonic International Flight Research (HIFiRE) Direct- 
Connect Rig (HDCR) - Mach 5.8 Flight Case 



O Expt, Body-Side Wall 

O Expt, Cowl-Side Wall 
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Hypersonic International Flight Research (HIFiRE) Direct- 
Connect Rig (HDCR) - Mach 8.0 Flight Case 
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Conclusions 


• BSL and SST models provide solutions on either side of experimental data 
for Shock-Wave / Turbulent Boundary Layer Interactions (SWTBLIs). 

• This work investigated alternative values for the shear stress limiter 
constant, a v which is set to 0.31 in the original SST model. 

• Incompressible models were investigated: For the Driver axisymmetric 
diffuser problem, SST provides best C p results; but perhaps fortuitous - 
considering turbulence measurements. 

• For SWTBLI problems, increasing a x results in less limiting of turbulent 

shear stress smaller separations. 

• Experimental data indicates values for a x = -u'v' / k larger than 0.31 are 
warranted - in agreement with computational results. 

• a x = 0.355 is recommended value for SWTBLI problems. 
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